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ABSTRACT 

When  non- conventional  guns  are  to  be  considered  or  when  detailed  design 
information  is  required,  Interior  ballistic  calculations  become  more  difficult 
and  time-consuming.  To  deal  with  these  problems,  the  equations  which  describe 
the  interior  ballistic  performance  of  guns  and  gun-like  weapons  have  been  pro¬ 
grammed  for  the  high-speed  digital  computers  available  at  the  Ballistic  Research 
Laboratories.  The  major  innovation  contained  in  the  equations  derived  in  this 
report  is  the  provision  for  use  of  propellant  charges  made  u|)  of  several  pro¬ 
pellants  of  different  chemical  compositions  and  different  granulations.  Re¬ 
sults  obtained  by  the  method  described  in  this  report  compare  favorably  with 
those  of  other  Interior  ballistic  systems.  In  addition,  considerably  more  de¬ 
tail  is  obtained  in  far  less  time.  A  comparison  with  experimental  data  from 
well- instrumented  gun-firings  is  also  presented  to  demonstrate  the  validity  of 
this  method  of  computation. 
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INTRODUCTION 


The  interior  ballistician  must  frequently  predict  the  interior  ballistic 

performance  of  guns.  In  some  instances,  it  is  sufficient  to  calculate  muzzle 

velocity  and  maximum  chamber  pressure  for  a  conventional  gun  from  a  knowledge 

of  the  propellant  charge,  the  projectile  wei^t,  and  the  gun  characteristics. 

(l)* 

This  calculation  is  usually  referred  to  as  the  classical  central  problem 
of  interior  ballistics.  When  non -conventional  guns  are  considered  or  when 
detailed  design  information  is  required,  it  is  necessary  to  know  more  than 
these  two  salient  values.  For  the  more  demanding  problems,  complete  interior 
ballistic  trajectories  may  have  to  be  calculated.  These  trajectories  consist 
of  displacement,  velocity,  and  acceleration  of  the  projectile  and  chamber  pres¬ 
sure,  all  as  functions  of  time. 

The  literature  of  interior  ballistics  contains  descriptions  of  many  methods 
for  solving  the  problem  of  predicting  the  performance  of  guns.  Methods, 

varying  from  the  purely  empirical  to  the  "exact"  theoretical,  have  been  devised 
in  tables,  graphs,  nomograms,  slide  rules,  and  simplified  equations  solved  in 
closed-form.  Some  of  these  methods  require  data  from  the  firing  of  the  gun 
being  considered  or  from  a  very  similar  gun.  All  of  these  methods  require  some 
simplification  of  the  basic  equations  of  interior  ballistics. 

To  eliminate  the  restrictions  imposed  by  assuirptions  made  only  to  facil¬ 
itate  the  mathematical  solution  of  the  problem,  the  interior  ballistic  equations 
have  been  programmed  for  hi^-speed  electronic  computers.  Both  analog  and  dig¬ 
ital  computers  have  been  used  to  calculate  detailed  interior  ballistic  trajec¬ 
tories.  There  are  advantages  and  disadvantages  associated  with  each  type  of 

(3) 

computer.  Several  years  ago,  the  interior  ballistic  equations  were  pro¬ 
grammed  for  the  digital  computers**  available  here  at  the  Ballistic  Research 
Laboratories.  Since  that  time,  considerable  use  has  been  made  of  this  program 
for  studying  gun  and  gun-like  systems  and  for  routine  calculations. 


*  Superscripts  Indicate  references  listed  at  the  end  of  this  report. 

**  Althou^  the  Interior  ballistic  equations  were  originally  programmed  only 

(4) 

for  the  ORDVAC,  they  have  been  recently  reprogrammed  in  more  general 
form  for  the  ORDVAC  and  the  newer  BRLESC. 
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The  computer  program  described  in  this  report  has  been  designed  to  solve 
a  set  of  non-linear,  ordinary  differential  and  algebraic  equations  which  simu¬ 
late  the  interior  ballistic  performance  of  a  gun.  In  this  method,  the  usual 
set  of  equations  which  pertains  to  the  burning  of  a  single  propellant  has  been 
modified  to  account  for  the  burning  of  composite  charges,  l.e.,  charges  made 

up  of  several  propellants  of  different  chemical  compositions  and  different 
* 

granulations.  The  computer  program  may  be  suitably  modified  to  study  non- 
conventional  guns  and  gun-like  systems.  A  number  of  these  optional  programs 
have  been  devised  and  used  extensively. 

INTERIOR  BALLISTIC  THEORY 

Interior  Balllatic  System 

The  basic  conponents  of  the  Interior  ballistic  system  for  a  conventional 
gun  are  shown  in  Figure  1.  A  set  of  equations  can  be  formulated  which  math¬ 
ematically  describes  the  distribution  of  energy  originating  from  the  burning 


Figure  1.  Basic  Conponents  of  the  Interior  Ballistic 
System  for  a  Conventional  Gun 


*  Wie  present  program  can  be  operated  with  as  many  as  five  different  types  of 
propellant  charges  for  each  problem. 

**  See  Section  entitled  Options  to  Routine. 
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propellant  and  the  subsequent  motion  of  the  components  of  the  system.  In  the 
development  which  follows,  two  major  assumptions  are  made  to  account  for  the 
behavior  of  composite  charges: 

1.  The  total  chemical  energy  available  is  the  simple  sum  of  the  chemical 
energies  of  the  individual  propellants. 

2.  The  total  gas  pressure  is  the  simple  sum  of  the  "partial  pressures" 
resulting  from  the  burning  of  the  individual  propellants. 

Energy  Equation 

Application  of  the  law  of  conservation  of  energy  leads  to  the  energy  equa¬ 
tion  of  interior  ballistics.  This  may  be  written  as: 

Energy  Released  _  Internal  Energy  ^  External  ^  Secondary 

by  Burning  Propellant  of  Propellant  Gases  Work  Done  Energy  Losses 

on  Projectile  (l) 


or: 


Q  =  U  +  W  +  Losses 


(la) 


In  Equation  (la)  the  energy  released  by  the  burning  propellant  (Q)  is  assumed 
to  be  equal  to  the  simple  sum  of  the  energies  released  by  the  individual  pro¬ 
pellants  as  previously  stated.  Therefore: 
n  T 


(2) 


Because  of  gas  expansion  and  external  work  performed  in  a  gun,  the  gas  temper¬ 
ature  is  less  than  the  adiabatic  flame  temperature  (T  ).  The  internal  energy 
of  the  gas  (u)  is  then: 


i=l 


The  external  work  done  on  the  projectile  is  given  by: 


W  =  A 


(3) 


(M 
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Substituting  Equations  (2),  (5),  and  (!»•)  into  Equation  (la)  gives: 

n  T  ”  /T  X 


which  may  be  rewritten  as: 
n  T 


dx  +  Losses 


(5) 


As  the  c  do  not  vary  greatly  over  the  temperature  ranges  from  T  to  T  , 
^i  °i 

they  can  be  replaced  with  mean  values  (c  ).  Integration  of  Equation  (5) 

,  ^i 

gives : 


n  X 

ZC.z.c  (T  -  T)  =  a  /  p,  dx  + 

•'o  '’ 


Losses 


(6) 


i=l 

and  solving  for  T: 

n 


T  = 


ZVi  VV  - 

i=l  ^ 


dx  -  Losses 


(7) 


U 

Zvi  V 

i=l  ^ 


Next,  the  force  of  each  propellant  is  defined  by: 


m.RT 
i  o. 


and  the  well-known  relations: 

SPl  -  5v^  .  m^R 

and:  c 


are  introduced. 


(8) 

(9) 


(10) 
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Combination  of  Equations  (9)  and  (lO)  gives: 


-1) 


Substitution  of  Equation  (ll)  into  Equation  (8)  gives: 


°i  (r  -  1)  c 

^  ""i  (12) 

Finally,  substitution  of  Equation  (12)  into  Equation  (7)  gives  Rfesal's  equation 
in  the  form; 


P.C.z 


^  — T - A  /  p,  dx  -  Losses 

Z.  ’'i  -  ^  ^  ^ 


P.C.z, 


V  1  1  i 


For  most  problems,  it  is  convenient  to  assume  the  igniter  completely  burned 
(z^  =  1)  at  zero- time.  Equation  (13)  may  be  restated  as: 


L  1=1  7^-1  -I  7i-l  0 

y  Vi!i  1  .  J 
(7,-l)To  J  ( 


p^  d  -  Losses 
b  X 


’'o. 


The  terms  A  ^  Losses  of  Equation  (l4)  can  now  be  considered 

in  more  detail.  The  work  done  on  the  projectile  results  in  an  equivalent  gain 
in  kinetic  energy  of  the  projectile  except  for  losses.  Including  these  losses 
under  the  general  category  of  energy  losses: 


-  X  W  p 

A  d*  .  1/2  -  V 
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According  to  Hunt,  '  '  the  energy  losses  to  be  considered  are; 

(1)  kinetic  energy  of  propellant  gas  and  unbumed  propellant, 

(2)  kinetic  energy  of  recoiling  parts  of  gun  and  carriage, 

(5)  heat  energy  lost  to  the  gun, 

(4)  strain  energy  of  the  gun, 

(5)  energy  lost  in  engraving  the  rotating  band  and  in  overcoming  friction 
down  the  bore, 

and 

(6)  rotational  energy  of  the  projectile. 

For  discussion  of  each  type  of  secondary  energy  loss,  see  Reference  (2). 

Types  (2),  (4),  and  (6)  are  estimated  to  be  less  than  one  percent  for  each 
category  and  have  been  neglected  here. 

The  kinetic  energy  of  propellant  gas  and  unburned  propellant  can  be  repre¬ 
sented  by 


The  energy  losses  resulting  from  heat  lost  to  the  gun  can  be  estimated  by  a 

(2) 

semi-empirical  relationship  described  by  Hunt: 


At  the  present  time,  the  introduction  of  a  more  sophisticated  treatment  of  heat 
loss,  with  its  attendant  complexity,  does  not  seem  to  be  warranted.  Such  a 
substitution  can  be  made  if  and  when  it  appears  desirable. 
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The  final  energy  losses  to  be  considered  here  consist  of  those  resulting 
from  engraving  of  the  rotating  band,  friction  between  the  moving  projectile 
and  the  gun  tube,  and  acceleration  of  air  ahead  of  the  projectile.  Individual 
estimates  of  these  are  difficult  to  make,  so  they  have  been  grouped  as  resis¬ 
tive  pressure  in  the  form: 

X 


(18) 


The  p  versus  x  function  is  discussed  in  greater  detail  in  the  section  concern- 
r 

ing  forces  acting  on  the  projectile. 


Substitution  of  Equations  (15)>  (l6),  (l?)^  and  (l8)  into  Equation  (l4) 
results  in  the  form  of  the  energy  equation  used  in  this  computer  program: 


r  ^  F  C,Z  T 

F^C^ 

2  /  V  P  \ 

+  ^  \  - 

^1  ^  ^i  ) 

Lit!  ’'i  J 

’'I 

2g  \W  +  i=l  / 

dx 


T  = 


F.C.z, 


Vi 


y  I  + 


(19) 


Equation  of  State 

The  pressure  acting  on  the  base  of  the  projectile  can  be  calculated  from 
a  series  of  equations,  once  the  temperature  of  the  gas  is  determined  from  the 
energy  equation.  Generally,  the  equation  of  state  for  an  ideal  gas  takes  the 
form; 


pji  = 


m^  RT 


(20) 


where  =  the  volume  per  unit  mass  of  i  th  propellant  gas. 

Now,  define  V^,  the  volume  behind  the  projectile  which  is  available  for  pro¬ 
pellant  gas,  as: 
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Volume  Available 
for  Propellant 

Gas 

Initial  Empty 

=  Chamber  Volume  + 

Volume  Resulting 
from  Projectile 
Motion 

Volume  Occupied 
-  by  Unburned 

Solid  Propellant 

Volume  Occupied 
by  Gas  Molecules 
( CO  volume) 

n  n 

or;  V  =  V 

c  0 

I  ■=!  (1--)  -S  VA 

i=l  ^  1=1 

__  / _ \ 

(22) 

By  the  definitions  of  Equations  (20)  and  (2l), 


^1  = 


(25) 


Substituting  Equations  (8)  and  (25)  into  Equation  (20)  and  rearranging  gives: 


P 


i 


F.C.z,  T 
111 


V  T 


c  o 


1 


(24) 


If  the  are  assumed  to  be  constants  over  the  tenrperature  range  from  T  to 
,  and  if  the  total  gas  pressure  is  taken  as  the  simple  sum  of  the  "partial 
pressures "  resulting  from  the  burning  of  the  individual  propellants  as  pre¬ 
viously  stated,  then; 


n 


T  “  Vi''! 


i=l 


c  i=l 


As  before,  if  it  is  assumed  that  the  igniter  is  completely  burned  (z^ 
at  zero-time,  Equation  (25)  may  be  restated  as; 

„  I-  n 


P  = 


I  '  o .  ’ 

1=1  i  I 


(25) 

1) 


(26) 
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The  space-mean  pressure,  p,  given  by  Equation  (26)  Is  used  in  the  cal¬ 
culation  of  the  fraction  of  propellant  burned  at  any  time.  This  relationship 
is  discussed  in  the  section  concerning  burning  rates.  There  is,  however,  a 
pressure  gradient  from  the  breech  of  the  gun  to  the  base  of  the  projectile 
which  must  be  considered  in  developing  the  equations  of  motion  for  the  pro¬ 
jectile.  This  pressure -gradient  problem  was  first  considered  by  Lagrange  and 
is  commonly  referred  to  as  the  Lagrange  Ballistic  Problem.  Later  studies  in 
this  area  were  made  by  Love  and  Pidduck,  Kent,  and  others.  For  this 
coirputer  program,  the  improved  Pidduck-Kent  solution  developed  by  Vinti  and 
Kravitz  has  been  used: 


Pb  = 


1  + 


n 

i=l 

W  S 
P 


(27)* 


In  addition  the  breech  pressure,  p^,  is  calculated  by  the  method  contained 
in  Reference  (6).  This  is  the  pressure  usually  measured  in  experimental  inte¬ 
rior  ballistic  studies: 


P„  = 


°  (i-%) 


-n‘-l 


(28) 


where :  l/ a 
'  o 


2  n'+3 
6 


2  (n*+l) 


i=l 


(28a) 


In  Reference  (6),  the  determination  of  6  depends  on  the 
heats,  y.  For  composite  charges,  an  effective  value 
is  used  for  this  purpose. 

7  '  = 


ratio  of  specific 
n 


i=l 


n 


Zc, 
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(27a) 


1 


and  n’  =  - 

7»-l  (28b) 

Mass -Fraction  Burning  Rate  Equation 

Both  the  energy  equation  (Equation  (19))  ajad  the  equation  of  state  (Equa¬ 
tion  (26))  are  algebraic  equations  whose  solutions  depend  upon  the  solutions 
of  several  non-linear,  ordinary  differential  equations.  The  mass-fraction  burn¬ 
ing  rate  equation  expresses  the  rate  of  consumption  of  solid  propellant  and 
hence  the  rate  of  evolution  of  propellant  gas.  This  may  be  written  as: 


For  certain  propellants,  including  those  plateau  and  mesa  types  used  in 
solid-fuel  rockets,  Equation  (52)  will  not  suffice  for  gun  calculations.  In 
these  cases,  it  is  preferable  to  make  use  of  a  tabular  listing  of  r. 's  and 
corresponding  p*s  (Equation  (50))  and  to  interpolate  for  the  desired  r^.  The 
r^'s  calculated  by  either  Equation  (50)  or  Equation  (52)  are  closed  chamber 
burning  rates.  As  discussed  in  later  sections  of  this  report,  these  burning 
rates  may  be  increased  by  addition  of  factors  proportional  to  the  velocity  and 
displacement  of  the  projectile  in  the  following  manner: 


r '  =  r  +  K  v  +  K  x 

1  i  V  X 


(52a) 
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Similarly,  the  form  function  described  by  Equation  (51)  may  be  stated 
in  one  of  several  ways.  In  many  Interior  ballistic  systems,  the  form  function 
is  chosen  for  convenience  of  analytical  solution.  Where  routine  numerical 
computations  are  handled  by  use  of  a  hi^-speed  digital  computer,  the  geomet¬ 
rical  form  of  the  propellant  grain  may  be  used  to  obtain  the  functional  rela¬ 
tionship,  f^,  between  and  z^.  For  the  usual  grain  shapes  encountered,  these 
equations  are  given  in  Appendix  A.  This  Appendix  also  contains  the  method  for 
handling  such  equations  in  the  computer  routine.  To  extend  these  equations  to 
Include  propellant  slivering  see  Reference  (9) • 

Equations  of  Projectile  Motion 

The  translational  motion  of  the  projectile  down  the  gun  tube  may  be  cal- 
cxilated  from  the  forces  acting  on  the  projectile.  Figure  2  shows  the  axial 
forces  considered  in  determining  the  resultant  force. 


F 

F 


f 
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Figure  2.  Axial  Forces  Acting  on  Projectile 

The  propulsive  force,  F^,  is  that  resulting  from  the  pressure  of  the 
propellant  gas  on  the  base  of  the  projectile  according  to: 

Fp  =  (55) 

where  p^  is  obtained  from  Equation  (27). 

The  frictional  force,  F^,  is  the  retarding  force  developed  by  resistauace 
between  the  bearing  surfaces  of  the  projectile  and  the  inside  of  the  gun  tube. 
This  is  usually  the  resistance  between  the  rotating  band  and  the  rifling  of 
the  tube  and  includes  the  force  required  to  engrave  the  rotating  band.  It  may 
be  expressed  as: 

=  P/  (5^+) 
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The  determination  of  Is  difficult  in  most  cases.  Many  interior 
ballistic  solutions  use  an  increased  projectile  mass  (approximately  5^)  to 
account  for  its  effect.  There  are  several  disadvantages  Inherent  in  such  a 
treatment.  Althou^  the  muzzle  velocity  may  be  calculated  reasonably  well, 
the  detailed  trajectory  will  be  altered  considerably.  It  is  not  possible  to 
simulate  the  case  where  a  projectile  lodges  in  the  bore  (see  Reference  (lO) 
for  experimental  trajectories  for  this  condition).  For  this  computer  program, 
experimental  data  of  the  type  given  in  Reference  (ll)  may  be  used  by  insert¬ 
ing  a  tabulation  of  the  function: 


P 


r 


G(x) 


(54a) 


The  gas  retardation  force,  P^,  is 
of  air  or  gas  ahead  of  the  prajectile. 


that  which  results  from  the  pressure 
stated  as: 


P  = 
r 


(55) 


where  p  is  small  enough  to  be  neglected  except  for  very  hi^  velocity  systems, 

S 

li^t  gas  guns,  and  other  special  applications.  In  the  discussion  of  the  Energy 

Equation  in  the  Interior  Ballistic  Theory  Section,  p  was  considered  a  part  of  p  . 

g  ^ 

The  resultant  force  in  the  axieil  direction  is  then: 


F  =  P 
a  p 


(56) 


or: 


F  = 
a 


(57) 


The  acceleration  of  the  projectile,  by  Newton's  second  law  of  motion, 


is 


or: 


A(p,  -  p  -  p  ) 
'^b  ^g  ^r' 


a  = 


M 

Ag  (p^  -  Pg  -  Pj.) 


a  = 


W 


(58) 


(59) 
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Summary  of  Interior  Ballistic  Equations 


The  equations  which  are  used  in  the  computer  progrsun  are  now  svumnarized 
for  ease  of  reference. 
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Equation  of  State 
r-  n 


-  T 
P  =  V 


(S^) 


L  1=1  “1 


T 

°I  J 


n  C. 


n 


where:  V=V+Ax-  Y'  —  (l-z  )  -V'  C  z  b 
c  o  /Lj  ^  1'  /Li  1  i  1 


1=1  ^1 


1=1 


^  ~  n 

1=1  ^ 


1+ 


W  S 
P 


= 


“  u-%) 


-n'-l 


Mass -Fraction  Burnlng-Rate  Equations 


(26) 


(22) 


(2?) 


(28) 


•  ^  J- 

dt“  (T  ^1  ^1 

% 

r,  =  P,  (P)  «1 


or: 

r’  =  r .  +  K  V  +  K  x 
1  1  V  X 


(29) 

(52) 


(52a) 


2K 


Equations  of  Projectile  Motion 


(59) 


(40) 


(41) 


COMPUTATION  ROUTINE 

The  set  of  non-linear,  ordinary  differential  and  algebraic  equations, 
summarized  at  the  end  of  the  previous  section,  simulates  the  interior  ballistic 
perfomiEuice  of  a  gun  or  gun-like  system.  A  numerical  computation  routine  has 
been  devised  for  the  simultaneous  solution  of  these  equations.  The  general¬ 
ized  flov-diagram  for  the  routine  is  presented  in  Appendix  B.  Using  the  FORAST 
(5) 

language,  '  '  the  solution  has  been  programmed  for  the  ORDVAC  and  BRLESC  com¬ 
puters. 


Prellmlnaiy  Routine 


To  reduce  computation  time  and  conserve  memory  space,  a  preliminary  routine 
has  been  introduced.  Here  all  data  required  for  the  computation  are  read  into 
the  computer,  constant  groupings  (e.g., 


F^Ci 


etc.,  are  calculated  and  stored 


for  subsequent  use,  and  data  to  permanently  identify  the  computer  run  are  print¬ 
ed  out.  A  complete  listing  of  required  input  data  may  be  found  in  Appendix  C. 
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Main  Routine 


The  main  computational  routine  is  presented  in  the  generalized  flow- 
diagram  of  Appendix  B.  To  follow  the  procedure,  consider  the  three  sequential 
phases  of  the  problem: 


Phase  I  -  Prom  time  of  ignition  until  the  projectile  starts  to  move. 

Phase  II  -  From  time  of  initial  projectile  motion  until  all  propellants 
are  consumed. 

Phase  III  -  From  time  of  propellant  burnout  until  projectile  leaves  the 
gun  muzzle. 

At  the  time  of  Ignition  (Phase  I  begins),  it  is  assumed  that  the  igniter 
is  completely  burned  =  l)  and  none  of  the  other  propellants  have  started 
to  burn  (all  z^  =  O).  The  space-mean  pressure,  consisting  only  of  the  igniter 
pressure,  is  calculated  from: 


P  =  Pt  =  Tr — 

^  \  (1^2) 
Equation  (1+2)  is  derived  from  Equations  (l9)  and  (26)  by  means  of  the  sin^)!!- 
fying  ignition  assuB^stions  stated  above. 


The  linear  burning  rate  for  each  propellant  can  now  be  determined  from 
either  Equation  (30)  or  Equation  (32)  in  combination  with  Equation  (32a).  If 
the  Interpolation  indicated  by  use  of  Equation  (30)  is  selected,  the  general¬ 
ized  interpolation  sub-routine*  is  employed.  The  mass-fractions  burned,  z^*s, 
during  a  small  time  interval,  dt,  are  determined  by  integration  of  Equation 
(29).  The  surface  areas  of  the  unburned  propellant  (see  Appendix  A)  are  used 

in  this  initial  calculation.  The  Runga-Kutta  method  of  numerical  integration, 

(12) 

as  modified  by  Gill,  is  commonly  used  for  the  solution  of  sets  of  or¬ 

dinary  differential  equations  and  has  been  employed  here. 

Calculation  of  the  temperature,  T,  from  Equation  (19)  and  the  volume 

available  for  propellant  gas,  V  ,  from  Equation  (22),  will  allow  the  calcu- 

c 

latlon  of  the  new  space-mean  pressure,  p,  at  time,  dt,  from  Equation  (26).  The 
surface  areas  of  the  now  partially  burned  propellants  are  computed  from  equa¬ 
tions  presented  in  Appendix  A.  All  results  of  Interest  are  printed-out  at  this 
time  **  and  these  results  used  as  initial  conditions  for  calculations  during 


*  See  Reference  (18)  for  interpolation  by  divided  differences. 

**  See  Appendix  C  for  listing  of  output  data. 
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the  ensuing  time-interval.  Those  terms  in  Equations  (ij);  (l9)j  and  (22)  which 
involve  velocity  or  displacement  are  zero  during  this  phase  of  the  computation. 
This  calculation-loop  is  continued  until  the  space-mean  pressure  exceeds  a 
pre-selected  "shot -start"  pressure  and  the  projectile  starts  to  move.  Phase  I, 
which  has  been  arbitrarily  defined,  ends  at  this  time. 

Phase  II  requires  the  addition  of  the  equations  of  motion  to  the  sequence 
followed  during  Phase  I.  Equations  (27),  (59)^  aad  (l^■l)  are  used  to  cal¬ 

culate  the  values  of  the  acceleration,  velocity,  and  displacement  of  the  pro¬ 
jectile  at  the  end  of  each  time  interval.  Integration  specified  in  Equations 
(4o)  and  (4l)  is  again  performed  by  the  Runga-Kutta-Gill  method.  Values  of 
velocity  and  displacement  are  now  available  for  use  in  terms  of  Equations  (l?)^ 
(19),  and  (22).  To  compute  values  for  Ep  =  A  ^  dx,  which  is  one  of 

the  terms  in  Equation  (l9)^  the  generalized  Interpolation  sub-routine  must  be 
used  to  obtain  p^  from  the  tabular  Information  described  by  Equation  (54a). 

This  integration  is  performed  by  use  of  the  Trapezoidal  Rule.* 

As  time  is  increased  by  the  addition  of  small  time -intervals,  calculations 
during  Phase  II  are  continued  around  this  expanded  loop  with  print -out  of  appro¬ 
priate  results  at  the  end  of  each  time  interval.  One  at  a  time,  the  propellants 
are  completely  consumed  and  this  phase  is  ended.  A  series  of  switches  has  been 
incorporated  in  the  program  to  circumvent  the  necessity  of  introducing  propel¬ 
lants  in  any  special  order.  In  fact,  it  may  not  always  be  possible  to  predict 
the  exact  order  in  which  several  different  propellants  will  be  burned  out.  The 
combination  of  the  propellant  switches  and  the  Start -of -mot ion  switch  makes  it 
possible  to  handle  problems  where  one  or  more  propellants  bum  out  before  the 
projectile  starts  to  move. 

With  all  propellants  consumed.  Phase  III  begins.  The  mass -fractions  burned 
have  all  become  unity  and  the  equations  concerned  with  burning  (Equations  (29), 
(31),  and  (50)  or  (32))  are  eliminated  from  the  loop.  As  in  the  other  phases. 


*  Althou^  the  Trapezoidal  Rule  is  a  relatively  crude  method  for  numerical 
integration,  the  accuracy  of  the  p^  versus  x  data  available  does  not  warrant 

a  more  acciirate  and  hence  more  complex  method. 
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results  are  printed-out  at  the  end  of  each  time -interval.  A  continual  check 
is  made  of  the  displacement  of  the  projectile  to  determine  whether  or  not  it 
has  reached  the  muzzle  of  the  gun.  When  the  projectile  passes  the  muzzle, 

Phase  III  has  ended  and  the  program  is  stopped. 

It  is  possible  for  the  projectile  to  reach  the  muzzle  (and  the  program 
stopped),  before  Phase  II  is  conipleted.  This  would  simulate  a  gun-firing  in 
which  unburned  propellant  is  ejected  from  the  muzzle.  It  is  adso  possible 
for  the  program  to  simulate  a  firing  in  which  the  projectile  becomes  lodged 
in  the  tube.  In  this  case,  Phase  III  is  not  completed  and  the  program  is 
stopped  when  the  projectile  displacement  does  not  increase. 

At  each  time-interval  after  the  beginning  of  Phase  II,  the  breech  pressure 
is  determined  from  Equation  (28)  and  printed  out.  This  result  is  not  used  in 
the  computational  routine  but  is  used  to  compare  theoretical  and  experimental 
results.  A  continual  check  is  made  of  the  calculated  pressures  and  the  maxlmiom 
breech  pressure  is  stored  with  its  associated  time  and  projectile  displacement. 
This  information  is  printed-out  at  the  end  of  the  program.  Calculations  dur¬ 
ing  the  last  time-interval  result  in  a  projectile  displacement  somewhat  greater 
than  the  desired  distance  to  the  muzzle.  A  linear  interpolation  between  re¬ 
sults  at  the  last  two  time -intervals  is  used  to  obtain  results  exactly  at  the 
muzzle.  These  results  are  also  printed-out  at  the  end  of  the  program. 

Options  to  Routine 

A  considerable  number  of  options  have  been  designed  and  coded  for  special 
problems.  These  Include  changes  which  enable  the  program  to  be  used  for  guns, 
or  gun-like  weapons,  which  are  not  of  conventional  design  (Figure  l)  and  changes 
which  vary  the  treatment  of  some  of  the  individual  parameters.  It  is  expected 
that  the  number  of  such  options  will  increase  as  the  program  is  used  for  a 
greater  number  and  variety  of  problems. 

Typical  options  for  non- conventional  guns  are  those  for  gun-boosted  rockets, 
traveling-charge  guns,  and  light-gas  guns  of  the  adiabatic  compressor  type. 
Examples  of  options  for  varied  treatment  of  individual  parameters  are  those  for 
adjusted  burning  rates  (previously  mentioned),  inhibited  propellant  surfaces, 
delayed  propellant  ignition,  variable  time-intervals,  constant  resistive  pres¬ 
sure,  and  resistive  pressure  as  a  function  of  base  pressure. 
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DISCUSSION 


No  attempt  has  been  made  here  to  present  a  new  and  different  Interior 
ballistic  theory.  The  objective  was  to  devise  a  convenient,  flexible  scheme 
for  performing  the  tedious  numerical  calculations  required  to  obtain  detailed 
interior  ballistic  trajectories.  The  selection  of  a  program  for  hi^-speed 
digital  computers  has  made  it  possible  to  eliminate  most  of  the  simplifications 
of  theory  required  to  facilitate  mathematical  solutions  by  other  methods. 

The  theory  presented  as  the  basis  for  the  conputer  routine  is  well-known 
and  has  only  been  modified  to  account  for  composite  charges.  There  are  several 
problems  present  in  all  Interior  ballistic  calculations  and  these  also  prove 
troublesome  here.  For  exanple,  useful  propellant  burning  rates  are  not  gen¬ 
erally  available.  It  is  known  that  burning  rates  obtained  from  experimental 
firings  in  closed  chambers  are  usually  low.  The  results  obtained  from  limited 
gun-firings  by  the  authors  indicate  gvin  burning  rates  may  be  twice  closed 

chamber  burning  rates  \mder  certain  conditions.  As  previously  mentioned,  op¬ 
tional  methods  of  adjusting  closed  chamber  burning  rates  have  been  provided 
for  in  this  program.  One  such  approach  is  to  consider  the  burning  rate  to  be 
a  function  of  the  projectile  velocity  (and  possibly  a  function  of  the  projec¬ 
tile  displacement)  in  addition  to  its  known  dependence  on  pressure.  This 
method  results  in  the  use  of  closed  chamber  burning  rates  when  the  gun  chamber 
is  practically  a  closed  chamber  (v  and  x  are  effectively  zero).  When  the  pro¬ 
jectile  is  moving  at  hi^er  velocities  and  is  further  down  tube,  reasonable 
Increases  in  burning  rates  are  obtained  and  used.  Other  equally  important  dif¬ 
ficulties  are  associated  with  the  determination  of  reasonable  values  for  resis¬ 
tive  pressure  and  shot-start  pressure. 

Considerable  versatility  has  been  built  into  the  program.  Instead  of 
stopping  the  computation  at  the  end  of  Phase  III,  a  new  problem  can  be  auto¬ 
matically  read  into  the  computer  and  solved.  This  multiple-case  feature  can 
be  employed  to  advantage  for  any  number  of  additional  problems  during  a  single 
computer  run. 

Typical  interior  ballistic  problems  were  used  to  compare  results  obtained 
from  this  computer  routine  with  results  from  other  interior  ballistic  schemes. 
(15),  (l4),  and  (l5)*  The  agreement  was  generally  very  good  when  the  other 
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schemes  vere  fairly  sophisticated.  In  addition,  detailed  interior  halllstic 
trajectories  are  produced  in  considerably  less  time  than  it  takes  to  calculate 
maximum  pressure  and  muzzle  velocity  by  other  systems.  A  typical  computer 
solution  for  a  conventional  gun  takes  only  10  seconds  If  magnetic  tape  output 
is  used  with  the  BRLESC. 

Results  from  computer  simulations  have  also  been  coii5)ared  to  experimental 
data  obtained  from  well -instrumented  gun  firings.  To  demonstrate  the  adequacy 
of  the  computer  routine,  data  from  a  typical  105mm  Howitzer  firing  were  pro¬ 
cessed  by  the  method  described  in  Reference  (ll).  In  Appendix  D  these  experi¬ 
mental  results  are  compared  with  the  predicted  results  obtained  from  a  simu¬ 
lation  of  this  firing. 


(?ouJi  3;^.  (5b  OXA^ 

PAUL  G.  BAER 

^JEROME  M.  FRANKLE 
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FORM  FUNCTION  EQUATIONS 


Geometrical  Equations 

1.  Initial  Volimie  of  a  Propellant  Grain 


V  =  f  (D.^  -  N.d/)L 
g.  4  '  i  1  i  ^ 


(A-1) 


where:  V  =  volume  of  an  unturned  propellant  grain,  in.' 


=  outside  diameter  of  grain,  in. 

Nj^  =  number  of  perforations,  dimensionless 
d^  =  diameter  of  perforation,  in. 

=  length  of  grain,  in. 

2.  'Volume  of  a  Partially  Burned  Propellant  Grain 

-  “i)  =  f  [<“i  -  “i  -  “i  '^1 "  "i>"]  'h  -  “1^ 


where:  z.  =  mass -fraction  of  i  th  propellant  burned  at  a 
given  time,  dimensionless 

u^  =  two  times  the  distance  each  surface  has  receded 
at  a  given  time,  in. 


5.  Initial  Surface  Area  of  a  Propellant  Grain 

2 

i 

=  n  iiLi.  t-w.a.jijj.j  +  - 

g. 


D.  -  N.d  ■ 

S  =  n  (D^  +  N.d.)  (L.)  + 

gj  Li  1  ^  2  J 


(A-5) 


where:  S  =  surface  area  of  an  \mbumed  propellant  grain,  in. 

h 


4.  Surface  Area  of  a  Partially  Burned  Propellant  Grain 


^1-  ” 


{K  +  ^(11^  »i)  ][  -  "J 


(A-4) 


where:  S^  =  surface  area  of  partially  burned  i  th  propellant 
grain  at  a  given  time,  in. 
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Equations  for  Hevrton-Raphson  Method*  for  Finding  Approximate  Values  of 
the  Real  Roots  of  a  Numerical  Equation 


1. 


2. 


Rearrange  Equation  (A-2)  to  set  f(u^)  =  0: 

=  f  L^{N^-1)  -  2(D^+I(^d^)  j 

-  [2L^(D^  .  N^d^)  +  (d/  - 

"  Ll  -  N,d/)}  -  (1-z^) 

Differentiate  Equation  (A-5)  with  respect  to  u^: 


f.  (u  )  =  - J  i3(Ni-l) 

du^  ^  t  ^ 

-  2  [l^(N^-I)  -  2  (D^  +  N^d^)]  u^ 

-  ^  -  v/)]} 


u. 

1 


2 


3.  The  value  of  the  root  of  Equation  (A-2)  is  then: 

f(Ui) 

Vl  =  *^1  ‘  f'(u^} 


(A-5) 


(A-6) 


(A-7) 


where:  improved  value  of  the  root,  where  the 

first  estimate  of  the  root  is  u^. 


Procedure 


1. 


For  each  propellant,  determine  by  integration  of  Equation  (29).  In  the 

initial  calculation  of  each  z^.  Equation  (A-3)  Is  used  to  compute  each 

S  (S.  =  S  when  u  =  O).  For  subsequent  calculations  of  each  z  ,  Equation 
1  1  1  1 

(A-4)  is  used  with  u^  determined  as  described  below. 


*  See  Reference  (l6)  for  a  discussion  of  this  method. 
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2.  The  z^’s  obtained  from  Equation  (29)  are  used  to  compute  the  u^*s  from 
Equation  (A-7)  and  then  the  new  S^’s  are  computed  from  Equation  (A-4). 

In  the  initial  calculation  of  u^,  the  first  estimate  of  its  value  is  zero. 

Equation  (A-7)  is  used  to  calculate  the  liiproved  value,  With 

as  the  estimate,  Equation  (A-7)  is  used  again  to  calculate  a  further-improved 
value,  This  procedure  is  continued  until  the  improvement  is  less  than 

10'^  inch. 
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APPENDIX  B 
Coirrputation  Routine 

1.  Generalized  Plow  Diagram 

2.  iOBAST  Listing 


l|.i 


INTERIOR  BALLISTICS  PROGRAM  FOR  GUNS 
GENERALIZED  FLOW  DIAGRAM 


Interior  Ballistics  Program  for  Guns 
FORAST  LISTING 


BLOCiPHl-Ptigni  ACl-AC5)Kl-K-/')y]-Y7»Ti-TS)')n-U5)St-S5)Rl-P5> 


•TS)')n-U5)St-S5)Rl  -P5> 


BLOC(01-U’7)Cl-C5)Pi-F5)GAl-GA?>)C0V1-CnV5)Tni-TO5)RHO1  -HH05 )  BE  T  1 -BE  T  5) 
PL0C(nCll-DCZ5)  ANl -AN1  1,4  J 


C0NTALPl-ALP5)l)l-U‘i)(DPl-0Pi>>Li-L5)NPj  -NP^;)  yCl-XC?n  )HCl-MC5> 
tNTf:H(A.«EAD)AN3)l?)» 


READ-EoRMAT(01)-<WP)AM)VC)AP)PE)nEl )PPmAX)* 

RE  aD-E()RmAT  (Oi  )  -  (  C  1  )E  nGAi  )T0I  )X _ 

RE4D-E()RMAT(01)-(DT)N1)KV)kX)D)EVP)«  SET(J=0>* 

PrAD-FORMAT(Ol)-(XCl.J)PRl,J)«  C0llNT(9n)IN(J)SnT0(Hl.l)» _ 

FMTER( INTEGER)N1)N)*SET<  J=0)*  I M ( NCP=-2*N ) * 

REaD-E0RPA7(01)-(  Cl,  J)F1.  J)GA1,  J)crivi  .  J)  TOI  .  J»WH01  .  J)* 


READ-F04MAT<01)-(BETl.J)ALPliJ)D1 ,J)nPl . J ) LI . J ) NP 1 . J ) » 
C0llNT<  ,N)IN(  j)G0T0(B?)SrSFT<  J  =  ri)X 


ENTER(A.PUNCHJANl)l)*  EMEP  (  A  .  PUNCH  )  AN8  9  )  I )  » 
ENTER<  A,P|JNCH)AN9)]  )» 


PUNCh-F0PHAT(0:->)-<l><WP>XM)VU)AP)0EL)PF>PPMAX)<A>* 
_ _  ENTER ( A . PUNCH  )  ANe9 ) 1 ) *  EN TEP ( A . PUNCH  )  AM  7  )  2  )  * 


P(JNCh-poPmaT(0.?)-<1>(CPFI  )GAI  )  ron<rRMTERA>* 
PUNCH-F0RHaT(0A)-<1>(C1»J)F1»JIGA1,J)COV1.J»TO1 .vDRHOl, J)<A>* 


CoiiNT<  .N)  IN(  J1«010(B4)*  SET(J  =  ii)*  EnTFR  (  A  .  PUNCH  )  An89  )  1  >  * 
ENTER« A.PUNCH)AN33)1)» 


PUnCR-F0RhaT(0S»-<1>(RET1, J) ALPl. J>ni  . J)nPi. J>L1. J»NP1.  J><A>* 
C0(|M  I  .  N  )  rN<  J  )r,OTO(B3  .  t  >*  SET  (  1=0)*  ENTER  <  A  .  PUNCH  )  AN«9  )  1  >  * 


ENrERi A.PUNCH)aN41)2)» 

PUNCh-F0RHAT(06)-<1>(XC1»J)PH1,J)<a>* _ _ 

COUN  I  < '20  )  IN(  J  )GOTO(  B3 .2J*  SET(J=0)»  F  NTER  <  A  .  PUNCH  )  AN89  )  1 )  * 

ENTER(  A.PUNCH)AN;?7)2)!t  SET(SnP  =  8  1  8.1) ,  IPxQ  )STUCK  =  B)  8. 5>)» _ 

punch-format < 07 )-<l>(DT)Nl)KV>KX)FVP)n)<A>« 

BClBF!«Cl/(aA|-l)«  ACl=HCl/T01*CCI»FI*CI/T01*  EVM=EVP»12« _ 

BCl,  J  =  F1,  j»Cl.J/(rjAl,  J-l»»  ACl.JsflCl  ,,)/Tol.  J*  CCl,  J  =  F1,  J*C1  .  J/ 

CONTTOi,  J* _ _ _ 

DC1.J=C1, j/RHOl.J*  ECl. J=C1/ J*COVl, J*  FCl . J=NP1. j-1 » 

GC).J»t.  1>J(NP1.J-1)-2(D1.J*NP1,J*DP1.J)« _ 

HCl  .  J  =  2*L1.  J(ni  i  J+NPl,  J*l)Pl,  J>+<D1,  J**2-NPli  J  +  OP:  ,  J**2)* 

_  ICt<J=Ll.J(Dl.J**2-NPl.J*DPl,J**2)* 


JC1  .-1*3  , 1416* ICI. J/4*C0UNt (.N)IN(J)GOtO(B4.1)*SET(J  =  0)* 

CT=OilTPl  =  ns 

fP'l  *C1,J*GA1  .“j  +  TPi*  CT  =  C1,j*CT*  COiiNT(,  N  )  I  N  (  J  )  GUfoTfr.  1 )  K 

GAP»TP1/CT*SET(J*n)!liGAF  =  GAP /(GAP-1)  HE  PsCT/wP* _ 

EP1 »1+£P /DEL »TPl  =  l/( GAP-1  )»TF2  =  l/( ( 24TP 1 *3 ) /DEL  * ( 2 1 TP  1 *2 ) /E P > * 
MCTD=RP*CT/D£L*TP4=0* AGH=AP*38fi . 4/WP* 


2*EXP(&AF*L0G( 1-TP2) )»  TP1=EXP(1.5*L0G(D>>*  TP 2»E XP ( 2 . 17 84LOG 
CONT(D))* 


TP40CI, J*T01, J  +  TP4*  C0UNT( ,N)  IN( J)GOTO( 85.2)*  SET(J»0)* 


^CT-29e>  )/(  (  1+0, 6*1 


MCI »<  . 38*12*TP1( XM  +  VO/AP) ( TP A/CT - 296 > ) / ( (1+0, 
CONTEVM**2 ) * 


LEAH(  /,M0S  .AT(Kl)*  CLEAH(7) 
pb=pr*pbr»xi»alp*intpr»o*  T*DT» 
CLEAR(5)N0S.AT(U1)*  XLi5T  =  0*  PRLSTeO* 
Y3,J=1.1*  C0UNT(5)1N(J)G0T0(B6.1)*  SETCJcO)* 


y3.J=0»  C0UNT( ,N)lN(J)G0T0(e6.2)*  SET(j=0)* 
IF-IMT(N*l)aOTn<e7.5)*  SET(SW3-Dr3.1 )* 


*2 
=  3 


0  0 '44'^ 
0045  ■ 
0  0  46“'] 
0047  I 


isa 


B7.6 

B7.7 

SET(sw4«Ri4»tiaoTO(Be)« 

SrT(SW5«R14  »»GnT0(R8)« 

B7.S 

SET(SH6*B14)»G0T0(8B)» 

.  0066 

B8 

TPl»li» 

0QP7 

B8.1 

TPl*DCl,  J  +  TPl*  COUNK ,N)  IN( J)GOTO(Be.l >*  SrT(J*n)» 
PTsCl*EI/(VO-TPl)»  SET(SW1=815)SHfl=P14,5>«  PMAXspTK 

0066 

0069 

B9 

ENTEW(R.K.G.  )DT)2,N)H9)Y1)K1)01>*  GOTOC  ,Sw1)K 

Tr(  Y5>  =  HP.OTO(R9. 1  )*SrT(SWll»Bl  0>.|sO)»r.OTO(DRl  )« 

B9.1 

Y3*l»K3«U<S1=O*Hl»O»u3=O»SET(Swll=0lO)jen)«QOTO(DR3.l)« 

17 

BIO 

IF(  Y4>  =  1  )G0T0(  RJ  0.]  )*SET.(sr11  =  B11  ).J8l  )«QOTO(  DRl  )* 

0073 

BlO.l 

Y4=lkK4a0SS2B0«R2a0X(j4sUSSET(SwllsB11)J>l>XGOTO(,SR3)S 

74 

Bll 

0075 

itaB 

Y5=l*Kt>ai,\XS3=0*R'>  =  U*y5=C*SET(SWll  =  Bl  ?)  J»2)*GnTOl  >SW4>* 

IF(  Y6>  =  l)GOTO(H3  2.1  >»SET(SV!11  =  813)  J=?)«GOTO(DBl  )* 

76 

0  077 

psfl 

Y6  =  l»,K6«0*S4  =  n»R4s(i*U'6  =  U»SET(SwllxBl.S>j3  3)*GOTO(  .Sw5)* 
ir(Y7>=l)GoTo(PlP.l>XSET<SWll=Pl4)J»4)»GOTo(nRl)» 

78 

0079  ' 

Y7  =  l*K7»o*S5=a*PP=ii*<^7  =  0*SET(SWll  =  8j  4)  j  =  4)*Gf)T0{  ,Sw6>* 

8n  ! 

SET( JBG )*TPls0» 

0081  i 

814.1 

TPl *DCl 1 j< 1-Y3  . J  ) +EC1 . J* Y3. J  +  TPl*COUNT<  «  N ) I N<  J ) GOTO ( 814 . 1 ) » 

0  08  2 

VC=Vt,*Ap»Xl-Tpl«SET  (  Ja0)*Tpl«BCl* 

0  00  3 

B14.2 

TPlrpCl ,  J*Y3> J  +  TP1*C0UN1 (,N)IN(J>G0T0(R14.2)« 

0084  I 

SET(J=0)*TP2=ACI* 

HIXEnD 

814.3 

TP?»ACl,J*y3# J+TP2%COuRT ( , N > I N ( J ) GOTO( Bl 4 . 3 > * 

SET(  Jsf!  )*TEMP  =  (  TPl-ALP  )/TP?«TPl  cCri* 

0  066  i 
0087  • 

B14.4  TP1»CCJ  ,J#y3,  J  +  TH1*C0UNM  ,N»IN(  J)a0T0(B14,<»)«  0068 

_ SET(  J  =  n)«PT«TEMP*TPl/VC» _ fiOTOt  .SW9)« _ Oud? 


B14.5  ENTER(R.K.GDI* 

DRl _ Rl,  J  =  BKT1  .,j♦£xP(  ALPI.  J*LOg(PT)  )»UO  =  U3  ,.]»H1  =  FC1  ,  J«H^  =  GCli  J« 


0090 

00V2 


JHH4.IC1 , J*HS  =  JC1,0(1-Y3,J)»H6  =  LI, J*H7bD1, J* 
H8  =  DPt>J»H9»NP3  .  J«HlfieDCl.  JKHlJaJCI  ,  J»60T  0(  G  A  M?  >  « 


0093 

0094 


DR3  Rl, J=Rl, J*KV*K?*«X4XI» 

K3.JcSl«J*RliJ/DCl<J» 


1094 

95 


DR3.1 

DR4 


GOTO< ,S411)« 

PB  =  PT/EPl»PflRspB/EP?« _ _ 

Kl=AGw«P8-PR)»K2  =  Y3!»XI=y2« 


1095 

0096 


0097 


DR  5 
DR9 

ENTEH( D.D.  INI  XI >PR1XC1)KRi )20>3I1)3 )* 

delx«xi-xlst«sumi=pr+prlst« 

0100  ! 
0090  1 

IntPRxI OELxiSUMl l/2*INTPR*XLSTsXI*PRLSTsPRK 

KEEQi 

ALP*<  McTd*K2**27  77?>.6)+AP*INTPR+MGL*K?»*?*  G0T0»B14.  S)« 

0100  j 

BIS 

B15.1 

iriPUPCiaoTOlBlS.l )»  SET(SR8»0R4>SN1  «B16)« 

PRR»PT*  PBsPTk 

0101  j 
010?  ' 

816 

^  IF ( ri>0)GOTO(8l7)»lF(Y3>»l) AND( Y4>  =  1 ) AND<  Y5>*1) AND( Y6>«1) 
CONTANDI  Y7>»1  )G0T0(B16.1UG0T0(B17)* 

n 

SET(STUCK«822)* 

Xf  =  Xi/12*V  =  Y1/12*  AF»«1/12*SET(  J»0  )1<ST=n» 

2103 

0104 

bczi.  J»K3  .J*Cl,J»"Cnul^T(  .t»l>  lN<J)QOTO(Bl  7.1)  *SET<  J=0)« 

STsSl ,  J*STsCOiiNT  <  ,N  )  I N  <  J  )  GOTO  (  817 . 2  )  *SE  T  (  J«  0  )  t 

0105 

01U6 

817.3 

,IFf PeR<PPHAXlGOTO(P17.3)»  ENtER ( A . PUNCH ) AN73 ) 1 ) *  G0T0(NEWRN)« 
IFIPMAX>PBR)G0T0(B38)»  PMAX*PBR*  XPHAXsXI*  TPMAXst* 

1106  ' 
107 

816 

816.1 

60T0(,S«P)» 

ENTER! A.PUNCH)AN1)1)» 

108 

109 

wm 

ENTER! A. PUNCH  )AN89)1>«  TM«T*100  0« 

PUNCm-FORHAT  !08  )-<!>!  TM)  XnP8R)PT  )PB)  V)  AF  )<A>* 

110  , 
111  ' 

818.3 

PUNCH-FOPHAT!09>-<1>!TM>XI)XF)TEMP)VC)PR)ST)<0>* 
PUNCH-f0RHAT!O10 )-<!>( TM)X I )Y3,J)DCZ1 .J)R1.J)S1#j)<6>« 

112  I 
111?  1 

j 

COUNT! ,N)IN!J)G0T0!B18.3)X  SET! J«0)X 

2112 

k6 


I 


IF(Xl>XH)GOT0(B«:l)!lXlLS1=XI*VLSTev*pB|  ST=PB* 

R<R.K.GL>JI  _ _  _  _  _ _ _ 


VmaX=I(Xm-XIlST)(V-V)ST)/(XI-XIlST>)*V|ST* 
PqMAx  =  ( < XM-XILST ) (Pfi-PBLST )/( X  I -XllST ) ) ♦[ 


TPMAX  =  TPMAX*1Q()0»  1117 

R(A.Pt/MCH)AA;H9)l)«  H:nTFR(  A.PllMf:H)AM31  )1  )% _  11 B 


PUNCH-F0RMAT(011)-^1>( VMAX)P'1AX)XpMAX)TPMAX)PBMAX><A>*  119 

fir)TO(  nf.^pn  ) 


F.NtER(  A.PUNCH)ANb9)l»*£iXfFR(A.PU(wCH)  AR97)l)*fNT£:R(  A.PUNCH)AN01)1)«  12o 
_ yMAX  =  0»PHMAX^KTPMAX  =  TPMAX  +  infln*  _ 111 


PUNCH-FOHMAT (01 1 )-<!>( VMAX)PMAX  »XPHAX)TPMAx )PBMAX )<A>*  1^2 

_ GQT0_(P1)*  _ _ 111. 


T1=HV-U0*T23HB+U0  0127 

FUr,7e&4(lJ0**3*Ml-li0**'2*M?-U0tH3*HA)-HS*  __ 


FPU=  .  78  54(  3*t)0**2*Hl-2*U0tH2-H3  )* 

U01  =UU-Fll/FPUlt  _ 


IF-ABS((IJ01-UO)<=.00001)QOTO(FF4)*MO=UOI»  0131 

n  I  0 ( GAM? ) * 


Sl  =  3.1416((H'4-ll0)(H7-U0  +  H9(Ha  +  IJ0))-t-.S*T1**?-.5*H9 
T*Tl**?>*Sl,  J=SI*M1  .  J  =  UO*fiOTO(nR3)* 


FORM( 10-1 0) 1-7) 

FQkM(l?-^-<y)  1-1)1  P-3-1  i,n-3)3-?)l  )  1  ?- 6 -0  1  3-3  )  12-6-0 -9  1 _ _ : 


f 0R«( lP-2-9) 12-7-9) 3-2) 12-1-7 ) 3  - 1 3 ) 1  P-4 -6 -p5 )  134 

FORM(  1P-P-9M2-7-9)  3-P)  1.-- 1"  7  )  3- 4  )  i  p- ?- 7  >  3- P  )  IP- 4-6  )  3- 4  )  1  ?-  0-6- 20  ) 


FORM( 12-9) 3-1 ) 1 2-6) 3-4) 12-6) 3-4) lP-6 ) 3-4 ) 1 p-l-7 ) 3-5 ) 1 2-1-3-P3 ) 

F0RH(3-2n)12-3-a)3-ll ) 12-5-7-3P) _ ! 

FORM (12-7 >3-5 >12-1-3 >3-4)12-9) 3- 1)12-9 >3-5)12-4-6) 3-5) 12-1-7- 17) 

FOFIM(  1  P-P-H)  3-1  )  1  2-3-6)  3 -P)  1 2- 6  -  1  n  )  IP- 6  - 1,) )  1  P- 6  -  i  n  )  1  ?- 5- 9  )  3^  J _ 


C0NT12-4-10-9) 

3 

_ F0RH(l2-2-6)3-1  >1 

?-3-8)3-2 

)12-P 

-8)3- 

2)1 2-4-8  )3-?)lP-5-l 0 

112-4-7)3-3) _ 

COtv|TlP-5-9-10) 

-8)3-2)12-1-7)3-3)12-4-9)3-] )l?-3-8 
FORM ( 12-5-6 ) 3-5) lP-6-8  >3-4 ) 12- 3-8 ) 3-? )1  P-?-7)3-3)l?-6-9-P4) 
(  B1  ) 


6 

IPROJ.  WT. 

barrel 

chamber 

105  MM  HOWITZER 
MORE  AREA  P-K 

RO 

765  t 

SS  PRESS  MAX  GUN  PRESSURE  1 

1 

1  charge 

FORCE 

gamma 

Ml  propellant 

COvULiJME  flame  T 

FMP 

/ 

DENSITY  / 

BETA  alpha  O.U.  GRAIN  OlA.  PERf  GR.IENGTH  NO.  PERp. 

. _ _  RESISTANCE 


PROj.  travel  pressure 

_ _ MiSCELL  ANE'OUS  _ 

or  NO.  PROP.  KV  KX  EST.  MU7.  VEL.  DIAMETER 

: _  P  GRLATER  than  desired  max  PRESSURE 


IMUZZLE  VEL.  MAX.  PHESSyRE  X  AT  PMAX  T  AT  PmaX  M|jZ  PRESSURE 
1 


1  PROJECTILE  Stopped 

33.  81.  153.  13.77  4600.  3.024  50000. 


60 .  ou _ 1  nOM , _ 

.6325  3670150.  1.264  31.08  2433.  .0567 

.5079-03  .8497 _ .  04/8 _ .  0194 _ .2453 _ 1. 

2.1356  3670150.  l.?64  31.08  2433.  .0567  26  i 

.5079-03  .8497 _ .  1344 _ .0142  .3127  7. _ 27  ! 


ru  CM  cmI 


APPENDIX  C 

Input  and  Output  Data 


1.  Input  Data 

2, '  Output  Data 

5.  Sairple  of  Output  Pormat 
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1.  INPUT  DATA 


Program 


Units 

Symbol 

Gun  Constants 

Wel^t  of  Projectile 

lb 

WP 

Length  of  Gun  Tube 

in. 

XM 

Empty  Volume  of  Chamber 

in.  5 

VO 

Cross-sectional  Area  of  Bore 

In.^ 

AP 

Shot -Start  Pressure 

psl 

PE 

Pidduck-Kent  Constant 

dimensionless 

DEL 

Resistive  Pressures 

psl 

PR1,J 

Travel  of  Projectile  Corresponding  to  each 
of  20  Resistive  Pressures 

in. 

XC1,J 

Diameter  of  Bore 

in. 

D 

Propellant  Physical  Constants 

Wel^ts  of  Propellants 

lb 

C1,J 

Weight  of  Igniter 

lb 

Cl 

Densities  of  Propellants 

Ib/in.^ 

RH01,J 

Outside  Diameter  of  Propellant  Grains 

in. 

Dl,  J 

Diameter  of  Propellant  Perforations 

in. 

DP1,J 

Length  of  Propellant  Grains 

in. 

LI,  J 

Number  of  Perforations  per  Grain 

dimensionless 

NP1,J 

Number  of  Propellants 

dimensionless 

Nl. 

Propellant  Thermodynamic  Constants* 

Forces  of  Propellants 

in. -lb/ lb 

F1,J 

Force  of  Igniter 

in. -Ib/lb 

FI 

*  See  Reference  (l?)  for  these  data. 
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Program 


Units 

Symbol 

Ratios  of  Specific  Heats  of  Propellants 

dimensionless 

GA1,J 

Ratio  of  Specific  Heats  of  Igniter 

dimensionless 

GAI 

Co volumes  of  Propellants 

in. ^/ib 

C0V1,J 

Adiabatic  Flame  Temperatures  of  Propellants 

°K 

T01,J 

Adiabatic  Flame  Temperature  of  Igniter 

°K 

TOI 

Burning  Rate  Coefficients 

in.  1 

psl« 

BET1,J 

Burning  Rate  Exponents  (a's) 

dimensionless 

ALP1,J 

Burning  Rate  Velocity  Coefficient 

in. 

sec  in. /sec 

KV 

Burning  Rate  Displacement  Coefficient 

in. 

sec-ln. 

KX 

Miscellaneous  Constants 

Time  Interval 

sec 

DT 

Estimated  Muzzle  Velocity 

ft  / sec 

EVP 

Maximum  Allowable  Breech  Pressure 

psi 

PPMAX 
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2.  OUTPUT  DATA 


Identifying  Data 

The  complete  list  of  input  data  is  printed  out  to  permanently  identify 


the  computation. 

Units 

Program 

Symbol 

Trajectory  Data 

Time 

milllsec 

TM 

Travel  of  Projectile 

in. 

XI 

Travel  of  Projectile 

ft 

XF 

Breech  Pressure 

psl 

PBR 

Space-mean  Pressure 

psi 

PT 

Base  Pressure 

psl 

PB 

Velocity  of  Projectile 

ft/ sec 

V 

Acceleration  of  Projectile 

ft/ sec^ 

AP 

Temperature  of  Propellant  Gas 

o 

K 

TEMP 

Volume  behind  Projectile  available 
for  Propellant  Gas 

in.  5 

VC 

Resistive  Pressure 

psi 

PR 

Total  Surface  Area  of  Propellants 

.  2 

in. 

ST 

Mass -fractions  of  Propellants  Burned 

dimensionless 

Y5,J 

Mass  Burning  Rates  of  Propellants 

lb/ sec 

DCZl, J 

Linear  Burning  Rates  of  Propellants 

in. /sec 

R1,J 

Surface  Areas  of  Propellants 

in.^ 

S1,J 
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OUTPUT  FORMAT 


PRfl.l.  WT. 

RaRREI 

chamreR 

G((N  PPfSSURF 

33.00000 

81.00000 

153,00000 

13.77000  3.02400 

4600  . 

50U00. 

CHARGE 

roRcE 

GAMMA 

Ml  propellant 

COvOLdME  FLAMF  TFMP 

DfNSITT 

.04290 

.6425(1 

1152000, 

3670150. 

1.2500 

1.2640 

2000. 
31.080  2433. 

.  056700 

igniter 

2.13560 

3670150, 

1.264(i 

31.000  2433. 

.  056700 

BETA 

.0005079 

ALPHA 

.0497 

O.D.  grain 
.(147b 

uiA.  peRf  qR. length 
.0194  .2453 

NO.  PERF. 

1. 

.0005079 

.8497 

.1344 

.0142  .3127 

7. 

PRO.I.  TRAVEL 

RESISTANCE 

PRESSiiRC 

j 

.000 

.100 

4500  . 
45011. 

i 

.200 

.350 

4500  . 
4500  , 

j 

.500 

1.000 

4500, 

4500. 

2.000 

3.50(1 

4500. 

4500. 

1 

4.000 

4.250 

2800  . 
2600. 

- 

4.500 

5.000 

2350. 

1900. 

5.250 

5.500 

1650. 
1400  , 

1 

1000. 

10  0  0  . 

30.000 

40.000 

1000  . 
1000. 

1 

50.000 

60.000 

1  0  0  0  . 
1000  . 

• 

1 

miscellaneous 

NO.  PROP, 

2. 

KV 

.0000000 

KX  EST.  MUZ. 

.0000000  1500. 

VEL.  diameter 
4 . 1340 

1 

( 

J 

i 

1 

f 

I 

j 

» 

. 

' 

■ 

. 

55 


TMmOOO  XX«.  000  ,  00UUT'Ntfr2054.?5  -WC'lnA.lA?  PR*  .  0  dT*  4  (1 18 . 93 

“nit.  1  ri  n  n  vTi  ..inn  V1t«.  n  n  i  c<  DC#ls  4  .  A  ii  1  ffis  .1  n  O  1  1  .  n  .X 


■nn^.iooo  «.*.n00  VH*.Q0a6  0C«»13.618  Rt»,10?  32-2357.10 


.  500 


,  5000 


.9000 


1.0000 


.0000  2202.89 


.0041  22.706 


103.975 
.170 


.170  2361.68 


4022.07 


.000 

1164.02  1164.02 

1164.0? 

.00 

.  1 

.  0  0  n 

103.917 

■  0 

4J23, 

.000  ,0096 


.000  1995.16  1995,16  1995.16  .00 

0  2308. 


36.496  .272  2369.10 


1.1000 


.000  .0138 


45.729  . 340  2374.28 


■"■See  list  of  Output  Data  for  Program  Symbols  and  Units. 
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1 .2u00 


.000 


.  0162 


51.055 


2377 .^2 


.3  79 


1.3000 

.000 

2919.97 

2919.97 

2919.97 

.00 

.  0  0  0  0 

_ 1.3000 _ 

_ dum _ 

.  oooo _ 

2346.09 

in^.pss _ 

_ .J] _ 

4035.09 _ 

1.3000  .000  .0445  39.548  .42?  1654.38 

1.3U  0  0  .  0  0  0  .  0  1  86  _ _ 56,911 _ .422 _ 2250 .-Zil 


1.400  0  .  000  3302.38  3302.36  3302. 3fl _ .jLQ _ ,  Q  B  Q  0 


■WilM 

.  nooo 

.u51  4 

V  w  u  f-  «  V>  y 

2355.58 

43.918 

103.112 

.469 

.  0 

1653,19 

4037.66 

1.4000 

.000 

.  0217 

63.344 

.  469 

2384.46 

i^HEEEHl 

3728.28 

102.953 

■m 

.  00  00 
4040.50 

■KlilM 

.  0590 
.  0250 

48.690 

70.407 

.520 
.  5?0 

1651.87 

2388.62 

1 .6000 

.000 

4202 . 09 

4202.09 

4202.09 

.00 

.  0000 

1.6000 

1 . 6  iJ  0  0 

.000 
.  0,0  0 

.  0  0  0  0 

.  0674 

2371.43 

53,898 

102.777 
.  576 

.  0 

1650.41 

4043.64 

1.6000 

.000 

.  0286 

78,156 

.  576 

2393.23 

1.7000 

1.7000 

.000 

.000 

4726  .'68 
.  0000 

4/28.60 

2377.98 

4728.68 

1 02.582 

.00 

.  0 

.  0000 
4047.11 

^HEEH 

.0767 
.  0326 

59.574 

86,655 

.637 

.637 

1648.79 

2398.32 

1.8000 

.  000, 

5384,73 

5312.49 

5169.10 

1.49 

8990.2031 

.  0000 
.  06  70 

2383.69 

65.753 

102.382 

.704 

4500.0 
1647 .00 

4050.93 

1.0000 

.000 

.  0371 

95.972 

.704 

2403.93 

1.9000 

1.9000 

.  004 
.  0  04 

HnlEEOI 

5959.75 

2388.72 

5798.89 

102.162 

■■om 

EeuEBHI 

17452.175 

4055.14 

,  0983 
.  0420 

72.457 

1Q6.156 

.777 

.777 

1645.03 

2410.11 

2.0000 

.  0  09 

6765.89 

6675.12 

6494.96 

4.94 

26604.618 

2.0000 

2.0000 

.0  09 
.009 

.  0007 
.1108 

2393.04 

79.730 

101.946 

.856 

4500.0 

1642.85 

4059.76 

2.0000 

.009 

.  0474 

117.296 

.856 

2416.91 

2.1000 

2.1000 

.  016 
.016 

7565.33 

.0014 

7463.03 

2396.65 

7262.38 
101 .739 

8.08 
4500 . 0 

37115.875 

4064.63 

2.1000 

2.1000 

.016 

.016 

.1245 
.  0534 

87,594 

129.452 

.942 

.942 

1640 , 46 
2424.37 

2.2000 

.  0  28 

8444 . 01 

8330.73 

6105,88 

12.30 

48449.250 

2.2000, 

2.2000 

■■ 

101.556 

1.035 

“500.0 

1637.63 

4070.30 

.026 

■  1 1  HUM 

1.035 

2432.55 
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4  ,,  R  HM  .<n:.!  T  T7C  Q _ Dn 

V6‘' _ 

.  3  a  n  n 

.  11 4  6 

9406 .19 

9280.1311 

9029. 5.V 

17.70 

60859.566 

^  .  3  0  U  0 
2.3000 

.  1  46 
.046 

.  0  0  3  9 
.1560 

24.il  .  74 

10  5.169 

101.410 

1.1.3‘! 

45on . n 

1  6  ,1 4 . 9  4 

40/6 .44 

2 . 30U0 

.  06  72 

1  5  /  .  0  5  2 

1.134 

2441 ,50 

2 . 4  iJ  0  0 
2.4U00 

■MM 

HHSnii 

HH 

10314.94 

2403.13 

1. 1)  0  .<  6 , 5  4 
101.320 

24 . 38 
4500.0 

74389,969 

4083 . 05 

2.4000 

2.4000 

.0  71 
.0  71 

.1740 
.  1)752 

11,4.896 

172.598 

1.242 

1  .  242 

1631.78 
2451 . 28 

2.5000 

.105 

11593.13 

11437 . 5.9 

11 1>8 .80 

32.46 

89066.992 

2.5000 

2.5000 

mogi 

^ni!i9 

mam 

2  4  0  3 . 7  0 

125,241 

1  0 .1  ,  3  0  4 
1.357 

4  5  0  0  ,  0 
1628 . 32 

4090.24 

2 . 5  0  0  0 

■m 

.0840 

189.356 

1.357 

2461  .92 

,2.6000 

2.6000 

wm 

12820 . 27 
.01  2  5 

12648 . 27 
2403.38 

123) 6.89 
101.383 

4  2 . 0  5 
4500,0 

1 04894.83 
4096. 04 

,l5(i 

.150 

.2150 
.  0936 

136.1/7 
20/. 338 

1.478 

1.478 

1624.55 

2473.48 

2.7000 

.20/  . 

14134.69 

13945, 06 

13568.68 

53.26 

121848.41 

2.7000 

2.7000 

.207 
.'2  0.7 

.0172 
.  2382 

241)2. 12 
14/. 660 

IQI  . 579 

1  .  6  0  7 

4  5  0  II  .  0 
1620.45 

4106.47 

2.7000 

.  2  0  7 

mmn 

226.532 

1.60  V 

2486 . 01 

2.8000 

2.8000 

,27B 
.  27fl 

15531.66 
.  0232 

15323,30 

2399.83 

149)9.72 

101.916 

66 . 19 
4500.0 

139066 .91 
4115.55 

2.8000 

2.8000 

.278 

.278 

.  2633 
.1155 

159.622 

246.895 

1.742 
1.74  2 

1  61* , Cl 
2499 . 55 

IHEHIIi 

.  366  . 

16322.35 

80.93 

158847 . 23 

2 .9000 
2.9000 

.366 
.  366 

.  0305 

.  2904 

2396 . 46 
171.973 

102.421 

1.882 

4 5 DO  .  0 

1611.19 

4l25.il 

2.9000 

.  366 

.  1280 

268.348 

1 .882 

2514.12 

3.0000 

•i.OOOO 

.473 

.*73 

18537.66 

.  0394 

18288.95 

2391.94 

1.7795.  3;< 
103.123 

97.58 
45(10.  0 

176638.52 

4135.75 

3.0000 

3.0000 

.^73 
.  *  73 

.  3196 
.1415 

184.593 
290 . 767 

2.0  27 
2.02  7 

1606.01 
2529 .75 

3.1000 

.6  0.1 

20119.30 

19849,  .37 

19313.63 

116.20 

199038.6* 

3.1000 

3.1000 

.601, 

.6  01 

.  0501 
.  3508 

2386 .2(1 

19/  .336 

104.050 

2.175 

4f>00  .  0 
1600.43 

4146 .6  7 

3.1OO0 

.801 

.1561 

313.984 

2.175 

2546 , 44 

3.2000 

3.2000 

.7  53 
.7  53 

21726.45 
.  0628 

21436.9.3 

2379.19 

2(18  58 .35 
105. 235 

136.62 

4500.0 

219793 . 76 
4158,66 

3.2000 

3.2000 

■KyH 

— 

■BOR 

210.040 
33/ -78  0 

2.323 

2.323 

1594.47 

2564.18 

.  3.3000 

■n 

23026.53 

22406 . 98 

1 59 . 46 

240601.50 

3.3000  ' 
3.30  0  0 

Him 

■REm 

■Bil^R 

Knm 

4171.09 

3.3000 

.931 

,1888 

361.890  * 

2.471 

2582.96 

■ 

58 


3.4000 

1.137 

.  0947 

2361.21 

1.08 . 51  0 

4500.0 

4184.12 

3.4000 

1  .1  37 

.4566 

1 581  .39 

3.4000 

1.137 

.2068 

386.006 

WBM 

2602.74 

3.5000  .1.373  26471.99  26116.83  2541  1.94  210.69  280976.61 

3.5000  1  ,373 _ .1164  /350.23 ^ lin.667  4500.0 _ 4197...72. 


3.5000  1.373  .4958  245.905  2.755  1574.29 

3.5000  1  .  373  .^?60  409.784 _ 2.755  2623.44 


■HI 

wm 

■HI 

RHHRII 

RHRIH 

wggm 

HH 

3.6000 

3.6000 

1.643 

1.643 

.13  6  9 
.  536  /  , 

2337.93 

256.417 

113.21  7 
2.886 

4500,0 

1  566 . 83 

4211.82 

3.6000 

1.643 

.2464 

432.062 

p.oe/' 

2644 .99 

HSnSIIIH 

1.948 

1.94a 

29280.32 

.1623 

28887.48 
2324 . 36 

28107.81 
116,1 02 

269.31 

4500.0 

317198.79 

4226.35 

3.7000 

3.7000 

1 .948 
1.948 

.  5792 
.  2678 

265.874 

454.872 

3.008 

3.0  08 

1559.05 

2667.31 

3 . 8  u  0  0 

2.29  0 

3  0  494.1  n 

'  3)1084.98 

29272.96 

3  0 1  .  0  6 

332854.33 

HWiliiM 

2 . 29 0 

2 . 2  9  (j 

.1908 

.6232 

2309.60 

274.108 

119.625 

3.117 

4500.0 

1550.97 

4241.23 

3.8000 

2.290 

.2902 

475.458 

3.117 

2690.26 

3.9000 
_  3. 9  0  on 

2.671 

2.671 

31549-.96 

.2226 

31126.67 

2293.75 

30266.56 

123.545 

334.20 
4500 . 0 

346472.95 

4256,38 

HEHEffiEHl 

2.671 
2.6  71 

.6683 

.3135 

280.985 

494.300 

3.212 

3.212 

1542.63 
2713 . 75 

4 . 0  0  0  0 

3 .0  92 

32430.96 

31995.85 

31132.28 

368.50 

357836.18 

4.0000 

4 . 0  00  0 

3.092  ' 

3 . 0  92 

.  2577 
../144 

2276.91 
286.41 7 

127.981 

3.293 

4500.0 

1534.07 

4271.70 

4.0000 

3.092 

.  3377 

511 . 1 26 

3.29  3 

2737,63 

Pilim 

3.555 

3.555 

33126.32 

.2963 

32661.88 

2259.23 

317*^9 .  a  n 
132.956 

403.76 

4353.8 

366805. 04 
4287.10 

3.555 

3.555 

•  /  6 1 
.  36?6 

29U.358 

525.726 

3.357 
3.3  57 

1625.33 

2761.77 

4.2000 

4.061 

33634.57 

33163.32 

32267 . 70 

440.69 

385840 . 09 

4.2000 

4.2000 

4.061 

4 . 0  6  1  ■ 

.  3384 
.8085 

2241 .04 
292.834 

138.493 

3.406 

2710.8 

1516.45 

4302.49 

4.  (i6l 

,  .3882 

536.002 

3.406 

2786 . 05 

4.3000 

4 . 3 IJ  0  0 

4.614 

4.614 

33945,66 

.3645 

33490.25 

2221.95 

32586.35 

144.626 

479.72 

2241.9 

403930.04 

4317.80 

4.3000 

4.3000 

4,614 

4,614 

Ksm 

1507.47 

2810.33 

4.4000 

5,213 

34  06  3.01) 

33605.99 

32696.97 

519.69 

412771.2(1 

■  4.4000 
4r4000 

5,213 

5,213 

.  4344 

.  9036 

2202.25 

293.326 

151.390 

3,452 

1688,0 

1498.44 

4332.93 

4.4000 

5,213 

.4406 

554,662 

3.462 

2634.49 

59 


yii-S-  MM  Hn^.I-T.Z.£.ia - Hn  TASi. 


K.OlSl 

33997 . 55 

.33541 .42 

32636.14 

560.46 

420082.51 

4. 5U00 

4. 50 no 

.  4884 
.95h9 

2182.08 

291 .497 

158.81 n 
3.4  52 

1097.9 

1 489.41 

4347.8? 

4.5U00 

nmginii 

.  46  /3 

559.430 

3.4  52 

2858.42 

4  .^6000 
4.6000 

6.558 

6.558 

33762.85 
.  5465 

33309 .87 
2161.5? 

3241 0.84 
166 . 9n4 

601.73 

658.5 

423888.9)) 

4362.40 

4.6000 

4.6000 

6.558 

6.558 

.  997  / 

.  4941 

288.442 

561.529 

1480.40 

2682.00 

4.7000 

7  .  ^  0  S 

.12683.80 

.32245.30 

31375.00 

642.69 

414921.61 

4.7000 

4.7000 

7.305 

7.305 

.6067 

1:  .  0  0  fi  0 

21 35,99 
.000 

176 . 078 
.  0  0 )' 

374.8 

.00 

2904.82 

4.7000 

7.305 

.52(|B 

556.704 

3 . 38  0 

2904 .82 

4.8000 

4.8000 

6  .  ].  (1  0 
8.100 

31546.35 

.6/5u 

31123.11 

2110.52 

3  0  283 . 1.  n 
185.947 

682.55 
29  0.7 

402  767 .92 
2926.94 

4.8000 

4.8000 

8 . 1  0  0 
8.100 

1.0000 
.  5469 

.000 

544.755 

.non 
3.28  3 

.00 

2926.94 

4.9000 

8.98k? 

3 u 38 2. 3  2 

29974 . 7  0 

29165.68 

720.95 

367359 . 05 

1^^ 

8.942 

8.942 

.7452 

1  ,  0  0  0  0 

2085.45 

.000 

196.490 

.non 

446 . 6 

.  no 

2946.17 

4.9000 

8.942 

.  5724 

531.861 

3.18? 

2948 ,17 

5.0000 

5.0000 

9.829 

9.629 

29211*63 

.8191 

28819.72 

2060.93 

28041.87 

207.686 

757.55 

883.7 

368342.90 

2966.51 

5.OO0O 

5.0000 

9.829 

9.829 

1  .  0000 
.  5973 

.  OuO 
518.326 

.  0  0  0 

3.  nen 

.00 

2968.51 

5.1000 

1  n.7  59 

28052.88 

27676.52 

26929. 53 

792.17 

348393.83 

5.1u00 

5.1000 

10.759 

10.759 

.  8966 

1  .1)000 

2037.15 

.000 

219.504 

.000 

1000.0 

.00 

2987 . 98 

5.1000 

10.759 

.6215 

504,441 

mm 

2987.98 

6.2000 

5.2000 

11.729 

11.729 

26920.30 

.9774 

26559.13 

2014.27 

25842.29 
231 .916 

825.32 

1000.0 

333785.59 

3006.61 

6.2000 

5.2000 

11 . 729 
11.729 

1 . 00 uo 

.6451 

.000 

490.446 

.000 

2.877 

.  0  0 
3006.61 

5.3000 

12.739 

25820 . 59 

25474.17 

24786.6? 

857.08 

319601.38 

5.3000 

5.31100 

■~12T7T5~ 

12.73.9 

1.0616 

1.0000 

1992.15 
.00  0 

244.894 

.00)' 

1000. 0 
.  0  0 

3024.42 

5.3000 

12.  739 

.66  79 

mm 

3024.42 

5.4000 

5.4000 

13.785 

13.785 

2476  0 . 24 

1 . 1488 

24428 . 04 
1970.86 

23768.73 
258.41 7 

887.50 

1000.0 

305924.79 

3041.46 

5.4000 

5.4000 

13.785 

13.785 

1 . 0  0  0  0 

.6901 

.000 

462.663 

5.5000 

14.868 

23743.17 

23424.62 

22792.39 

916.62 

292806.55 

5.5000 

5.5000 

14.868 

14.868 

1.239U 

1 . 0  0  fl  0 

1950.39 

.000 

272.461 

.000 

1000.0 

.00 

3057.75 

5.5000 

14.868 

2.590 

3057.75 

6o 


1  u  t;  MM  hH  ]  ]  T  -/FP 


t?ll  /6!^ 


5. 6110  0 

IS.  98.6 

227/1 .6',! 

22466 .00 

21859.72 

94^.52 

28  0275. 0  3 

J, .  6  U  11  0 

5 . 6  mi  0 

is.'^Bs 

16.905 

1.332" 

’  .  'iomi 

I93U.73 

.non 

287. one 

.  '■  ri  !■ 

imiD .  n 
.  0  0 

3073.33 

6 . 6 11  0  0 

16.905 

.  7326 

435.883 

2.501 

3  0  7  3 . 3  3 

5. 7UQ0 

5 . 7  Cl  0  11 

17 .134 
17.134 

2) 6“6 . 4l 
1.4278 

21553.31 

1  9 1 1 . 8  7 

209  /1  .  58 
302.027 

971 .23 
1000.0 

266341.81 

3088.25 

5.7uno 

6.7000 

17.134 

17.134 

1  .  II  I)  It  Cl 
.  7529 

.000 

423.041 

.  0  0  0 
2.456 

.  0  0 

3  n  8  8 . 2  5 

6.0000 

1, 0 . 3 1  6 

2ij96  7.52 

20686.21 

20127.89 

996.83 

257005.85 
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